Submitted to Phys Rev E, July 4 1996 



Electroconvection in a Suspended Fluid Film: 
A Linear Stability Analysis 

Zahir A. Daya and Stephen W. Morris 

Department of Physics and Erindale College 

On ' University of Toronto 

0\ 



in 



X 



60 St. George St., 
Toronto, Ontario, Canada M5S 1A7 

John R. de Bruyn 



I Department of Physics 



Memorial University of Newfoundland 



o 

^ . 

Q> , St. John's, Newfoundland, Canada A IB 3X7 

O ; (February 9, 2008) 



Abstract 



A suspended fluid film with two free surfaces convects when a sufficiently 
large voltage is applied across it. We present a linear stability analysis for 
this system. The forces driving convection are due to the interaction of the 
applied electric field with space charge which develops near the free surfaces. 
Our analysis is similar to that for the two-dimensional Benard problem, but 
with important differences due to coupling between the charge distribution 
and the field. We find the neutral stability boundary of a dimensionless 
control parameter 7?. as a function of the dimensionless wave number k. TZ, 
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which is proportional to the square of the apphed voltage, is analogous to 
the Rayleigh number. The critical values Tic and found from the 

minimum of the stability boundary, and its curvature at the minimum gives 
the correlation length ^o- The characteristic time scale tq, which depends 
on a second dimensionless parameter "P, analogous to the Prandtl number, is 
determined from the linear growth rate near onset, and tq are coefficients 
in the Ginzburg-Landau amplitude equation which describes the flow pattern 
near onset in this system. We compare our results to recent experiments. 
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I. INTRODUCTION 



The regular patterns which form when a dissipative, nonequihbrium system is driven just 
beyond the threshold of certain symmetry-breaking instabilities are in many ways analogous 
to the simple ordered phases which appear following equilibrium phase transitions Pat- 
terns can, however, exhibit interesting nonlinear dynamical behavior, for example chaotic 
motion, which has no analog in equilibrium systems. Several fluid-dynamical systems un- 
dergo pattern-forming instabilities which are amenable to both theoretical and laboratory 
study. Examples which have been extensively explored into the nonlinear regime include 
Rayleigh-Benard convection ||l]J^, Taylor vortex flow and electroconvection in nematic 
liquid crystals In each of these cases, an essential foundation for understanding the 
nonlinear behavior is a complete analysis of the initial linear instability. The linear stability 
analyses for Rayleigh-Benard convection and Taylor vortex flow are classic problems in fluid 
mechanics The mechanism of the linear instability for electroconvection in nematic 

liquid crystals required many years to elucidate, but even this very complex system is now 
reasonably well understood in both the linear and the nonlinear regimes 0. 

Our objective in this paper is to carry out a realistic linear stability analysis for a different 
electrically driven instability, namely electroconvection in a thin, suspended fluid fllm. We 
have observed electroconvection patterns in experiments on thin suspended fllms of smectic 
A liquid crystals i,0,||,|,|lOl , which are isotropic in the plane of the fllm but have a layered 
structure which very strongly impedes flows perpendicular to the fllm. As a result, these fllms 
can convect rapidly with no change in thickness. We have observed convection in fllms only 
a few molecules thick. Our immediate goal is to understand the onset of electroconvection in 
this system. The instability mechanism we describe is, however, not speciflc to smectic fllms 
and would presumably apply to any sufficiently two-dimensional, weakly conducting fluid 
fllm. In fact, similar convective flows have been observed in thicker fllms of nematic and 
isotropic liquids |ll|] . In these cases, however, surface tension effects and the convective flow 
itself cause thickness variations in the fllms, which make their behavior more complicated 
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than that of the smectic films. 

The model we describe below is physically similar to a highly simplified one proposed by 



Faetti, Fronzoni and RoUa [|TT| for the "vortex mode" convection they observed in nematic 
films, but our analysis is carried much farther. There is also some similarity between the 
driving mechanism considered here and that which drives electroconvection near the free 



surface in a partly-filled capacitor |12 



The relevant experimental arrangement is shown schematically in Fig. |l[ A thin fluid 
film is suspended between electrodes, with both its top and bottom surfaces free. The width 
of the film d is much larger than its thickness s, and we will treat it as being purely two- 
dimensional. When the dc voltage applied across the electrodes exceeds a critical value, 
the film convects in a pattern of vortices confined to the plane of the film. We neglect any 
effects of air drag by assuming that the film is suspended in vacuum. We will also treat the 
film as a weak ohmic conductor, and neglect any electrochemical charge production on the 
electrodes or in the bulk of the fluid. The currents involved are assumed to be sufficiently 
small that magnetic forces are insignificant. 

The body force responsible for driving any electroconvective flow results from an electric 
field acting on regions of nonzero charge density in the fluid. To analyze the electroconvection 
system we must first identify the mechanism which gives rise to regions of locally unbalanced 
charge in the fluid, and second, solve for the charges and fields selfconsistently, since these 
are coupled by Maxwell's equations. In our model, the charge density arises due to the 
electrical boundary conditions at the free surfaces of the film The inset to Fig. |l] shows 
the essential details of the charge separation mechanism. Below the onset of convection, 
the applied voltage drives a uniform, steady current density J through the film. This is 
accompanied by a constant electric field Ej^gj^^g = J/cr, where a is the bulk conductivity. 
The interior field Ej^^gj^-^g has no component perpendicular to the film plane. However, the 
exterior field E^^^-gj^^g must have both parallel and perpendicular components just outside 
the free surfaces of the film. It cannot in general be perpendicular to the surface because 
the surface is not an equipotential: the film is an ohmic conductor, so its surface potential 



varies linearly with the coordinate between the electrodes. The parallel component E , .-, 
^ ±-1- outside 

is equal to Ej^^gj^^g by the usual matching conditions on electric fields across surfaces. The 
perpendicular component E^^^gj^jg is proportional to the surface charge density a at that 
location. It is the interaction of the parallel component of the field Ej^^gj^-^g = Ej^^^^^gj^g with 
the surface charge density a at the two free surfaces that drives the convective fiow above 
the onset of the instability. 

In Section we calculate the surface charge density below onset by solving for the 
fields exterior to the film. This problem is solved analytically for thin films in two simple 
electrode geometries. We find that a "charge inversion" is set up in the base state: the film 
has a positive charge density close to the positive electrode, and a negative charge density 
close to the negative electrode. This inverted charge distribution is sustained by the apphed 
potential difference across the conductor — without a potential difference, the film surfaces 
are equipotentials, the component of the field parallel to the surface is zero, and thus there 
are no forces to drive convection, even if an electrostatic surface charge is present. This 
inverted base state configuration is analogous to the mass density inversion that arises in 
the Benard problem, and in Section |IT1| we show that it leads to a hydrodynamic instablity 
when the applied voltage is sufficiently large. Unlike the density in the Benard problem, 
however, the charge density in the base state is a nonlinear function of position across the 
film. 

By treating the film as two-dimensional, we neglect the diffusion of charge on the scale 
of the film thickness s, which acts to smear the surface charge over a thickness of order 
the Debye screening length [§]. For the very thin films considered here, \d may be 
comparable to s, in which case the surface charges and surface forces described above will 
extend over the whole thickness of the film. One can show that the total charge contained 
in one such Debye layer is the same as that which would reside at one surface in the absence 
of diffusion. The approximation that the film is a two-dimensional conducting sheet may be 
expected to break down for thick films for which s ^ A/5. In this limit, surface forces may 
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lead to significant shears and internal flows, as is apparently the case in thick nematic and 
isotropic films [[T^]. Diffusion on the much larger scale of the film width d is also neglected. 

In Section |T|, we describe the linear stability analysis for infinitesimal perturbations 
about the base state. The stability calculation is somewhat analogous to that for the Benard 
problem. Dimensionless quantities TZ and V appear, which are analogous to the Rayleigh 
and Prandtl numbers. TZ is proportional to the square of the applied voltage, while V is the 
ratio of the thin film charge relaxation time to the viscous relaxation time. 

The differences between our calculation and the Benard problem are due to the additional 
coupling between the field and the charge density, which is also responsible for the nonlin- 
earity of the base state charge density. The charge density and the field are analogous to the 
mass density and the gravitational acceleration in the Benard problem; the new requirement 
that these also satisfy Maxwell's equations amounts to requiring a nonlocal relation between 
the charge density and the electric potential. If we suppress this nonlocality by assuming 
these are simply proportional, the base state charge density becomes linear and our problem 
reduces completely to the Benard case. Interestingly, this proportionality is nearly correct 
except near the edges of the film. 

In Section |IV|, we discuss the results for the neutral curve, and compare the predictions 
for the critical voltage Vc and the critical wave number to the values obtained from exper- 
iments. We also calculate the correlation length ,^0 from the curvature of the neutral curve 
near k^. and the characteristic time tq from the linear growth rate at /tc. These quantities, 
which are coefficients in the Ginzburg-Landau equation which describes the amplitude of the 
pattern near onset, as discussed below are also compared to experimental results. 



Section is a brief summary and conclusion. 

II. THE BASE STATE CHARGE DENSITY 

Our first task is to calculate the configuration of charges and fields below the onset of 
convection. As described in the previous section, this is essentially an electrostatic problem 
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in the region exterior to the film. The coordinates used and the geometry of the electrodes 
and film are shown in Fig. ^ The origin is at the centre of the film, which lies between 
z = ±s/2 and y = ±d/2. We will consider the limit of a thin film for which s -C The 
film is assumed to extend infinitely in the x direction. The upper and lower surfaces of the 
film are free and the region outside the film has permittivity eg. The permittivity e of the 
fluid will turn out to be irrelevant to the analysis of the base state. 

The film will be treated as a charged conducting sheet of negligible thickness in the xy 
plane with bulk conductivity a. Its edges at y = ±d/2 are held at applied potentials ±V/2 
by electrodes of zero thickness. Below the onset of convection, the film behaves as an ohmic 
conductor, so that the potential on the film varies linearly between —V/2 and +V/2 for 
—d/2 < y < d/2. The potential is zero on the x axis and as \z\ oo. The potential exterior 
to the film is symmetric above and below the xy plane and independent of x. The charge 
density q is proportional to the perpendicular component of the field exterior to the film, 
at the film's surface, and so to the z derivative of the potential there. To calculate q we 
need only solve for the potential in the upper half of the yz plane, subject to boundary 
conditions on the y axis. We will consider two simple electrode geometries, which we refer 
to as "plates" and "wires". In the plate geometry, we specify the potential on the rest of 
the y axis to be —V/2 for y < —d/2 and +V/2 for y > d/2. Solving for the potentials is a 
Dirichlet problem which we solve below using a Green function. This geometry corresponds 
to a film held between infinite knife edges. Most of the experiments []6|,|8|,PJT0[] , however, 
used thin wire electrodes to support the edges of the film, as shown schematically in Fig. 
|l|(a). To model this geometry, the applied potential ±V/2 is specified only at the two points 
y = ±d/2. For \y\ > d/2 we require that the z derivative of the potential on the y axis be 
zero. Thus, in the upper half of the yz plane, we must solve a mixed boundary value problem 
with Dirichlet conditions for —d/2 < y < d/2 and Neumann conditions for \y\ > d/2. This 



is done analytically below, using the method of dual integral equations P 
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A. The Base State for Plate Electrodes 



We begin with the simpler plate electrode configuration. We must solve the Laplace 
equation for the potential \l/ in the upper half yz plane, 

+ l^b(l/'^) = 0> (2-1) 



subject to the piecewise linear Dirichlet boundary conditions on the y axis 

^(z/,0) = -- - oo <!/<-- 

V d d , , 
- jy -2<y<2 (2.2) 

V d 

= — - < -u < oo. 

2 2 - ^ 

The appropriate Green function is constructed from a unit line charge at {y', z') and its 
image at {y', —z'), 

G{y,z;y ,z ) = -log- . 2.3 

[y - yr + + z'y 

The potential at any point in the upper half plane is given by 

=ir*fe'.o)|i #'=ir (2.4) 

AtiJ-oo oz' TxJ-ooz^ + [y - y ) 

The surface charge density on the film is a consequence the fact that the z components of 
the electric fields inside and outside the conducting film are different. Inside the conductor, 
in the absence of diffusion, the z component of the field is identically zero, as in Fig. [^(b), 
and hence only the external field is required to determine the surface charge density. As 
mentioned in the Introduction, the presence of a diffusion layer near the surface does not 
change the total charge density present, per unit area. The surface charge density on upper 
side of the film is —eod'^ /dz\z=s/2- In the limit s — > 0, d"^ /dz is discontinuous across z = 0, 
so we use a one-sided derivative valid as z ^ 0^. To get the total surface charge density q 
on the film, we introduce a factor of two to account for the two free surfaces, so that 



oz 



(2.5) 



Using this with Eq. (|2.4|) , interchanging the order of differentiation and integration, and 
using Eq. (p.2|) gives 



TT 

'2^ 



-V/2 



+ {y- y') 



+ 



Vy'/d 



+ 



=0+ •'-2^'^ + i.y-y') 



=0+ 



(2.6) 



=0+ 



+ {y- y'f 

where the subscript p denotes plate electrodes. After integration, the resulting expression 
was expanded as a power series in z and evaluated as 2; ^ 0^. After some simplification, we 
find the surface charge density for the case of plate electrodes to be given by 



2^60 



log 



y-d/2 



y + d/2 



(2.7) 



This distribution is shown in Fig. ^ Note that the charge density is positive near the positive 
electrode and negative near the negative electrode, giving the charge inversion described in 
the Introduction. The charge density diverges at the electrodes, which is an unphysical 
consequence of the limit s 0. In the real system, the finite thicknesses of the film and 
electrodes will impose a cutoff on qp. This divergence, while unphysical, is weak enough to 
be mathematically tractable. It will turn out that the rigid boundary conditions we impose 
on the flow, described in the next section, ensure that the contributions from the edges of 
the film are small. On the other hand, the fact that qp{y) is not a linear function of y has 
important consequences for the quantitative results of the stability analysis. 

B. The Base State for Wire Electrodes 



We now turn to the case of wire electrodes. For this case, the mixed boundary value 
electrostatic problem can be solved by the theory of dual integral equations ||14|. We must 
solve Eq. ( p.l|) , subject to the mixed boundary conditions 



V 

^{y,0) = -y 



d^{y,z) 



dz 







=0+ 



for \y\ < d/2 
for \y\ > d/2. 



(2.8) 
(2.9) 
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By separation of variables and using the fact that \I/(0, 0) = 0, we make the ansatz that the 
potential in the upper half plane can be written as 

^(y, z) = r ^^e-^' sin iky)dk. (2.10) 
Jo k 

With this ansatz, we find the dual integral equations 

r sin (ky)dk = ^y for \y\ < d/2 (2.11) 
Jo k d 

/ A{k) sin {ky)dk = for \y\>d/2. (2.12) 
Jo 

This pair of integral equations may be solved for the potential in the upper half plane by 
the method of Sneddon 



giving 



V rMkd/2)^_,, 



^(y, ^) = Y X \ ^ sin {ky)dk, (2.13) 

where Ji is the first order Bessel function of the first kind. Using Eq. ([^.5|), the total surface 
charge density is given by 

qM = eoV / Mkd/ 2) sin {ky)dk = , ^ • (2.14) 

•^0 d ^{d/2f - y^ 

The subscript w denotes wire electrodes. This result is also shown in Fig. |^. As for the 
plate electrodes, we find an inverted charge distribution and divergences at the edges of the 
film. 

Below the onset of convection, the electric field inside the film is constant and points 
along — y. The force driving the convection is due to the in-plane electric field acting on the 
charge densities q. Unlike the Benard problem, in which the temperature profile is linear 
below onset, here the charge density, and hence the body force, is not linear in y. This 
has the effect of introducing certain non-constant coefficients into the stability problem, as 
described in the next section. 



III. LINEAR STABILITY ANALYSIS 



In this Section we consider the stability of the base state to infinitesimal perturbations. 
We will show that for sufficiently large V, the electrical forces overcome viscous and con- 
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duction losses and the film becomes unstable to convection. 



A. The Perturbation Equations 



Within the thin film, we assume that the fiuid velocity u is confined to the film plane, 
with u — u^c + vy. In addition, we will neglect any shears in the z direction. As discussed 
above, these assumptions are reasonable in the context of thin smectic films, where the 
layer structure strongly inhibits fiow across layers. We treat the film as a two-dimensional 
conducting fluid, with areal material parameters Ps — sp, rjs — sr], as — sa, where s is the 
fllm thickness and p is the bulk density, rj is the bulk molecular viscosity, and a is the bulk 
conductivity. (In smectic A films, the viscosity is highly anisotropic; the relevant component 
to use for 77 is 773, the viscosity related to shears within layer planes.) The two-dimensional 
pressure field is given by Pg = sP. The fiow is driven by the surface force density gE^, where 
Eg is the electric field in the film plane. The fiow velocity is governed by the two-dimensional 
Navier-Stokes equation 



where Vs is the two-dimensional gradient xd/dx + yd/dy. The fluid is assumed to be 
incompressible, so that 



This condition may also be viewed as a constant thickness assumption. We impose physically 
realistic rigid boundary conditions on u at the edges of the film, so 




(3.1) 



V, • u = 0. 



(3.2) 



u = and ^ = at y = ±d/2. 



(3.3) 



The motion of charge is governed by the charge continuity equation 




dq 



(3.4) 
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in which is the two-dimensional current density in the plane of the film, and includes 
contributions from both conduction (o"sEs) and convection (gu). Diffusion of charge in the 
plane of the film has been neglected. 

The electric field in the plane of the film Eg is given by 

E, = -V.^s = -Vs^L^o • (3-5) 

As in the previous section, "if is the potential which solves the three-dimensional Laplace 
equation, 

V^^ = 0, (3.6) 

in the half space z > 0, with the surface charge density given by Eq. (|2.5| ). 

The surface charge density q in Eqs. (|3.1| ) and ( |3.4| ) in principle contains both the density 
of free charges q^^ee ^^^^ dielectric polarization charges, so that q = g£j.gg — V^-Ps. 

Here = sP, where P, the bulk polarization density, is given by P = cqxE, where x is the 
electric susceptibility. Inside the film, E = E^, independent of z, and has zero z component. 
Eq. ( |2.5| ), which only involves exterior fields, therefore only gives q^Yee' ■'■^ base state, 
however, • E^ = and q = q^^ee^ polarization effects are irrelevant. In the general 
case, one can show using the scalings given below that the dimensionless form of the ■ P^ 
terms are proportional to x^/d- For the experiments of interest, x^/d ^ 10"^, so we can 
safely neglect polarization effects. 

Eqs. (|3.1|) -(|3l6|) have a simple solution when ii = 0. Eq. (|3.6|) , subject to the appropriate 
boundary conditions, corresponds to the base state charge density problem solved in the 
previous section. Once the fields and q are found from Eqs. ( |2.5|) and (p.6|), Eq. ( p.l|) can be 
solved for the pressure Pg whose gradient balances the surface force density gE^. Eq. ( p.4|) 
then gives a constant current density = c^Eg, with a constant interior field E^ = —(y/d)y. 
This is a selfconsistent solution which corresponds to hydrostatic equilibrium. The pressure 
gradient points everywhere toward the midline of the film, ?/ = 0. 

Expanding the divergence in Eq. ( |3.4| ) and using Eq. (^.2[), the charge continuity equa- 
tion becomes 
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dq ^ ^ ^ 

— = -u ■ VsQ - osVs ■ E5 

ot 



(3.7) 



Note that Vs ■ is not equivalent to V ■ E, because of discontinuities in the z component 
of E at z = 0. Thus, Vs ■ E^ is not directly related to the charge density q. Vs ■ E^ can be 
found only after solving the full three-dimensional Laplace problem given by Eq. 
To examine the stabihty of the base state, we introduce perturbed quantities 



u = + uW, 
Es = Elo)+E«. 



(3.8) 
(3.9) 
(3.10) 
(3.11) 



where E^ = Ey^^^y and E« = E^^^^x + Ey'^'^^y. Here Ey'^^^ = -V/d, and is the base 
state charge density found in the previous section. To first order in the small perturbations, 
Eqs. (13), (13), and (|3) become 



Vs ■ u^^) = 



Ps 



dt 



-VsPs^^) + r^.Vs^fl^^) + q^'^Ej'^± + q^'^Ey^'^y + q^'^Ey^'^y 
-u(i) ■ Vsg - cxs [Vs ■ (^.('^x + E,(i)y) 



(3.12) 
(3.13) 

(3.14) 



Taking the curl of Eq. ( ^.13| ) eliminates the pressure. Taking a second curl and using Eq. 
( P^ gives 

d 



PsT^Vs'iJ^') = r/sVsVs'i3(') - Vs X 
ot 



Vs X (gWE.(i)x + g(i)i?,(°)y + q^'^Ey^'^y) 



. (3.15) 



From this equation, we select the y component, which is 

ot (7X^ Oy OX 

Using the fact that g'-^^ is only a function of ?/, Eq. ( p.l4| ) becomes 



dt dy dx dy 



(3.16) 



(3.17) 



13 



In the previous equations, E^^^ is the in-plane electric field produced by the charge distribu- 
tion q^^\ Introducing perturbations for the electric potential with = \E'(°) + \I'(^\ we have 
E^^-* = — Vs\E's^^'' = —Vs'^^^^^^Q- The perturbation potential \E'*-^'' is the solution of a new 
three-dimensional Laplace problem in z > analogous to Eq. ( p.6|) , 



V^*^^) = 0, (3.18) 

with = -2eo(9^(^V^^)l.=o+- 

We now replace the various field components with the appropriate derivatives of the 
potential in Eqs. ( p.l6|) and (|3.17| ), to get 



4VM^^ - MS'^^ + E^-A^ + '-l^'-lt (3.19) 



dt dy ^\ dx"^ dy"^ 

Eqs. (|3.18|) - (|3.2CI|) are the equations for the perturbations that we must solve to determine 
the stability of the base state. 

The specification of the boundary conditions necessary to solve Eq. ( p.l8| ) requires some 
explanation. By writing \E' = \E'*^°) + \E'*^^\ we split the full Laplace problem of Eq. ( p.6|) 
into separate Laplace problems at each order. At zeroth order, the boundary conditions at 
z = on \E'(°) are those described in Section || for each electrode configuration in the base 
state. In particular, \I''^'') was set equal to ±V/2 at the edges of the film. At first order, the 
boundary conditions on \E'*^^) require that "^^^^ = at the edges of the film, and on both of 
the electrodes in the plate case. In the wire electrode case we require d'^^^y dz\z=o+ = for 
\y\ > d/2 and z = 0. In either electrode case, we will find Dirichlet boundary conditions 
for '^^^^ on the film itself by selfconsistently solving Eq. ( p.20| ). Proceeding in this way, the 
boundary conditions on the total potential are satisified by the superposition of "i!^^^ and 
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B. The Normal Mode Expansion 



We now expand the velocity, charge density and potential perturbations in normal modes 
which are periodic in x with wavenumber k, and have growth rate 7, 



(3.21) 
(3.22) 
(3.23) 



We substitute Eqs. (|3.21| )- (|3.23| ) into Eqs. ( |3.18|) -( pl20D and non-dimensionalize the system 
by dividing lengths by d, times by t^djos and charge densities by e^^Vjd. We then write 
D = d/dy and define new dimensionless quantities k = kd and Q{y) = d'^Dq^^\y)/eQV. The 
resulting dimensionless equations are 



and 



(z}2 _ ^2) ( ^2 _ ^2 _ 2 ) A + n^nie - QQs 



{D^ - K'^)n, - 76 - gA = 0, 







where Qs = ^|z=o- The three-dimensional Laplace equation, Eq. (|3.18| ), becomes 



52 92 
+ 



Qy2 Q^2 

with Eq. ( p.5|) imposing the condition that 



oz 



n = o, 



(3.24) 



(3.25) 



(3.26) 



(3.27) 



z=0+ 



Eq. ( |3.26| ) is a two-dimensional Helmholtz equation in the half plane x = 0, z > 0, which 
is perpendicular to the plane of the film. Eqs. ( p.26| ) and ( |3.27|) determine the rather com- 
plicated nonlocal coupling between the in-plane potential function Qs{y, 1^, 7) = ^{y, 0, k, 7) 
and the charge density function Q{y, ^,7). 

Two dimensionless groups appear: TZ, which plays the part of the Rayleigh number, and 
V, which plays the part of the Prandtl number. In terms of the bulk, rather than surface, 
material parameters, they are given by 
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7^ = -^-^ (3.28) 

and 

V = . (3.29) 

pads 

TZ, the control parameter, is proportional to V"^. It is interesting to note that TZ is inde- 
pendent of d, the width of the film. The Prandtl-like parameter V may be regarded as the 
ratio Tg/r^ of the two time scales in the problem, the charge relaxation time for thin films 
@ T(j = eod/as and the viscous relaxation time = pd'^/rj. 

The effect of the nonlinear y dependence of the derivative of the base state charge density 
Dq^^\y) is contained in the non-constant coefficient Q{y)- For the plate and wire electrode 
configurations, we find from Eqs. ( |2.7|) and ( p.l4| ) that Q{y) is given by 



and 



= (1 _ 4^2)3/2 ' (3-31) 



respectively. 



C. Analogy to the Benard Problem 

The above equations bear a strong analogy to the corresponding equations in the Benard 
problem. The correspondence becomes complete if the nonlocal coupling of the charges and 
potentials given by Eq. ( |3.26| ) and ( |3.27| ), is suppressed by simply putting q oc "^g- Applying 



this assumption to the base state removes the nonlinear y dependence of the charge density 
so that Q{y) = 1. In fact, detailed analysis shows that q is always nearly proportional to 
in the central part of the film. This can be seen, for example, in Fig. ^ near y = 0. If this 
proportionality is assumed to hold everywhere, then our continuity equation for charge, Eq. 
( p.4|) becomes identical to the thermal diffusion equation in the Benard problem. Under the 
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same assumption, the force term gE,, in Eq. ( ^.11 ) becomes proportional to qy, which is the 
form of the analgous term in the Benard problem. Turning the argument around, one can 
say that the reason that our system does not reduce to the Benard problem is because the 
charges and fields are nonlocally coupled via the the charge distribution's own self-field. 



D. The Compatibility Condition 

To find the conditions which TZ, V, k and 7 must satisfy for solutions to exist, we solve 
the linearized equations by means of various expansions. A crucial step that must be done 
numerically is the solution of the Helmholtz equation in the plane perpendicular to the film, 
Eq. ( |3.26| ), which necessarily involves a numerical relaxation calculation. 

At the edges of the film, y = ±1/2, the rigid boundary conditions on the flow velocity 
u^^\ given by Eq. require that A(i/) satisfy the four conditions 



A(±l/2) = L)A(±l/2) = 0. (3.32) 
To ensure this, we expand A{y) as 

A{y) = A^Cmiy), (3.33) 



m=l 



where the Cm{y) are even Chandrasekhar functions 



cosh(A^y) cos(A^y) 
Cm{y) = — 1/^ • ^^-^^^ 

cosh(Am/2) cos(Am/2) 



Here A.^ is the mth root of tanh(Am/2) +tan(Am/2) = |T^ . We can restrict the expansion 
to even functions because of the symmetry of the equations about y = 0. Note that an 
expansion in Cm has been shown to give a good description of the velocity field measured 
in experiments on electroconvection in smectic films |^ . Only relative amplitudes matter in 
Eq. (|3.33|) , so we set Ai = 1. It follows from linearity that we can also write Q = Y,m Am^m 



and 6 = J^mAmQm, where Qm and 6^ are the solutions corresponding to A = Cm- As 
above, we denote flgm = ^mL=o- 
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1. The potential function ^sm for 7 = 



Setting 7 = and substituting Cm for A in Eq. (|3.25|) gives 



{D^ - K^)n,m = QCm, (3.35) 
which may be solved directly by Fourier expansion. Since QCm is even, we expand 

oo 

QCm = J2 ^rnn COS (2n7r?/), (3.36) 

n=0 

in which 

bmo = 2 r Q{y)Cm{y)dy (3.37) 
Jo 

and 



bmn = 4: Q{y)Cm{y) cos {2miy)dy . (3.38) 
Jo 

Using a similar Fourier expansion of Qgm and imposing the zero boundary conditions at 
y = ±1/2, we find 



E 



^0 [(2n7r)2 + 



-!)"■ cosh (ny) 



cos {2miy) 



(3.39) 



cosh (k/2) 

To calculate Qsm, we used a Romberg numerical integration scheme [16| to tabulate the 
integrals for bmn in Eqs. ( p.37|) and (|3.38|) for the each of the two electrode geometries, using 
Qp and Q^ as given by Eqs. ( p.30| ) and ( |3.31|) . We used an upper cutoff of n = 29, which 
was dictated by the double precision accuracy of the Romberg scheme. 



2. The potential function for 7 7^ 
For nonzero 7, we solved 

by an iterative scheme. We used the 7 = solution, Eq. ( |3.39| ), to find a first approximation 
f2sm'°'- From this, we calculated the corresponding approximate the charge density function 
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Qm using the relaxation algorithm described in Section [III D 3| below. Then QCm + 'jQ 



was Fourier expanded in the same manner as QCm in Eqs. ( p.36D -( p.38| ) above. This expan 



sion was used to find a series solution analogous to Eq. (p.39|) for the next approximation 
^sm^^\ which was then relaxed to find &m^^^- This sequence of steps was iterated until it 
converged for both fi^m aud G^- The convergence criterion was a sum of the squares of 
100 differences in successive iterates distributed on < ?/ < 1/2. For I7I < 0.1, the sum 
converged after 7 or 8 iterations to a precision limited by the Romberg integration scheme 
used to find the Fourier coefficients. 

3. The charge density function Qm 



We solved the Helmholtz equation, Eq. ( p.26|) , for Qm for each of the two electrode 



geometries, using a simple SOR algorithm |[T^. In each case, the Dirichlet conditions on 
Qm for —l/2<y< 1/2 and z = are given by Qsm (Eq. ( p.39| )) in the case of 7 = 0, or 
by the corresponding expression during iteration for 7 7^ 0. Beyond the film, for \y\ > 1/2, 
z = 0, we applied the Dirichlet condition = in the plate electrode case, and Neumann 
conditions {dQm/dz)\z=o = in the wire case. 

Because Qm is even in it need only be relaxed in the first quadrant of the yz plane. 
We used an x square lattice of cells in this quadrant, with NfUm < N points between 
y = and y = 1/2. On the outer edges of the lattice, we set Qm = to enforce the zero 
boundary condition at infinity. Starting with = 100 and NfUm = 50, we systematically 
increased and Njum is such a way that NfUm/N 0. All the quantities calculated below 
showed a small residual monotonic variation with NfUm', we removed this by plotting each 
against l/Nfn^ and extrapolating to l/NfUm — > 0. 

From the resulting Qrn, the charge density perturbation Bm was determined from Eq. 
( p.27| ) by taking the one-sided z derivative numerically. was therefore only known at 
NfUm lattice points across the film. For the purposes of integration, we used a Chebyshev 
interpolation [jl6l of these points. 
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4- The compatibility condition 



E 

m=l 



= 0. 



(3.41) 



To find the general compatibility conditions on solutions, we substitute the expansions 
for A, fig and 6 into Eq. (g]^) to get 

Multiplying by the Chandrasekhar function Ci{y) and integrating from y = —1/2 to 
y = +1/2, we form inner products, denoted by (• ■ •). Then Eq. ( p.41|) becomes a lin- 
ear homogenous system with the determinant compatibility condition 

K^n(^Ci[em-Q^s,n)) 

After some simplification, this can be written as 

2^ ' 0, (3.43) 
where Fim = {CiiQm — Q^sm))- The matrix elements Xim are given analytically by 0] 



0. 



(3.42) 



Xlm — {Ci Cm) 

2 



A; - A^ 



when / 7^ m 
when / = m. 



(3.44) 
(3.45) 

(3.46) 



where = D'^Cm{y), etc. The matrix elements F;m(/t, 7) were evaluated numerically for 
each electrode configuration using Romberg integration |16|. The divergences in Q{y) at the 
edges of the film are overcome because Ci{y) goes to zero sufficiently fast at y = ±1/2. The 
functions and Qsm are simple smooth functions for the first few values of m and are 
straightforward to integrate numerically. 



E. Marginal Stability 



To find the conditions for marginal stability, we set the growth rate of the perturbations 
7 equal to zero in the compatibility condition, Eq. ( p.43| ). The Prandtl-like dimensionless 
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group V drops out, so that the marginal stabihty conditions are independent of P, just as 
is the case in the Benard problem. Eq. ( |3.43| ) then implicitly defines the marginal stability 
curve TZ = 71o{k). We proceeded as follows. Choosing a value of k, we set / = m = 1 
and calculate Fii{k). Then Eq. ( |3.43| ) can be simply solved to get the first approximation 
7?.W(k). We then find Fim{K) for l,m = 1,2 and search near 7?.W(k) for roots of the 2x2 
determinant, Eq. (|]4|), to find 7^[,2](/^). We can then use = 1 to find A2 in Eq. (p.33|) 
by backsubstitution. We carried this algorithm to the third order in the Chandrasekhar 
expansion, for which the maximum value of 1 is of order 10~^ and the resulting neutral 
curve TZo{n) no longer changes significantly. Fig. | shows the amplitudes A2 and ^3 for 
the wire case, relative to = 1. It is clear that the higher terms in the Chandrasekhar 
expansion contribute very little to the sum in Eq. ( |3.33|) . 

Fig. ^ shows the neutral curve for the plate and wire electrode configurations. The 
minima of these curves define the critical values Kc and TZc = TIo{k,c) for each case. These 
values are listed in Table |. We find that both neutral curves give Kc between 4 and 5, but 
that TZc is lower for the wire electrode case. This is apparently due to the steeper slope of 
q^'^\y), evident in Fig. ^ for the case of wire electrodes. Neither value of is particularly 
close to the Benard value of 3.117, but they are in reasonable agreement with the value 
determined from the smectic film experiments [|6t^|lOi, as discussed in section below. 

We can define a length scale ^0 in terms of the curvature of 7lo{n) near llJl)i]: 



^° " 2 dK^ 



(3.47) 



where ec = {TZo{k,) /TZc) — 1- This length will appear as a coefficient in an amplitude equation 
description of the convection pattern near onset ||l|,0,|l. To find ^0 accurately, we fit Cc to 
a parabola over a range k = ± Ak and then systematically reduced An until the value 
of .^0 taken from the fit became independent of An. This corresponded to a fitting range 
< 5 X 10^^ The values of ^o, 

given in Table |I|, were slightly dependent on the electrode 

configuration. 
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F. The Linear Growth Rate 7 



Returning to the full compatibility condition Eq. (|3.43|) with 7 7^ 0, we consider the 
behavior of the growth rate 7 of the linear modes near the critical values of TZ and k. The 
time scale tq, defined by 

.-1 57(e) 



(3.48) 



where e = {TZ/TZc) ~ 1, will also appear in an amplitude equation description of the pattern 
near onset |]T],0,§|- 

The matrix element Fimin, 7) is rather expensive to calculate for 7 7^ 0, because we must 
use the iteration scheme outlined in Section [111 D 2| . It is most computationally efficient to 



choose a value of 7, fix /t = Kc, and then solve Eq. (|3.43| ) for TZ. This was done for ten 
values of 7 in the range —0.1 < 7 < +0.1, using three Chandrasekhar modes. The results 
depend on V. The resulting function 7(e) is very nearly linear in e with a "P-dependent slope 
and 7(0) = 0. We determined tq from polynomial fits to 7(e) for V > 0.01. The results are 
only slightly dependent on electrode configuration, tq is plotted as a function of V for wire 
electrodes in Fig. 0. For V > 1, tq tends towards the limiting values given in Table |. 



IV. DISCUSSION 

The stability analysis presented above demonstrates that a thin, weakly conducting 
suspended fiuid film becomes unstable to spatially periodic convective fiow if a sufficiently 
large voltage is applied across the film. Since our analysis is linear, it cannot describe the 
convection pattern above onset, but it does provide important information about the onset 
of convection. In this Section, we discuss the theoretical results in the light of previous 
experiments on smectic [P,[7|,p|,p!0| and nematic |jll| films. 

Eq. ( p.28| ), combined with the neutral curve, predicts that the onset of convection occurs 
at a critical voltage proportional to the film thickness s, and independent of the film width 
d, given by 
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V, = -Jar]n,. (4.1) 
Co 

The dependence of Vc on y^or/ follows inevitably from dimensional analysis. Faetti et al. also 
found a critical voltage proportional to s^Jorj and independent of d from a highly simplified 
model of the "vortex mode" observed by them in nematic films |jlT . 



In experiments on convection in smectic films, Vc has been found to be proportional 
to s for films up to about 20 molecular layers (i.e., about 63 nm) thick For larger 

s, Vc grows somewhat more slowly. This may be a sign that layer-over-layer shears in the 
z direction exist for thicker films; such flows are not accounted for in our calculation. A 
linear dependence of Vc on s has also been observed in experiments on nematic films []TT|. 
The nematic films are much thicker than the smectic films, and have significant thickness 
nonuniformities. The also exhibit slow flows even below the onset of convection, making 
Vc{s) rather difficult to measure. 



Our most recent experiments on smectic films [10| show no dependence of Vc on d for 
films with d between 0.7 and 2.0 mm, with thicknesses s between two and 25 molecular 
layers, that is, between 6.3 nm and 80 nm. This is consistent with the prediction of Eq. 



(O). Over about the same range of thickness, as noted above, Vc is also proportional to s, 
as predicted. A weak variation of Vc with d was, however, observed in our earlier work [Q for 
d in the larger range of 0.36 mm to 3.5 mm. This work used a thicker film (107 molecular 
layers, or 340 nm) and a slightly different electrode configuration, with guard electrodes 
outside the main electrode wires. These features may have contributed some d-dependent 
three-dimensional effects. 

The wavenumber at onset observed in smectic film experiments is [|lOl fcj:^^* = 4.94 ± 
0.03 d~^ . The measured value of d is uncertain to ±5%, so this result yields a measured 
dimensionless wavenumber k^^^* = 4.94 ± 0.25. This is in satisfactory agreement with 
the value of Kc = 4.74 found from the minimum in the calculated neutral curve for wire 
electrodes. At present, no data is available for comparison to the predictions for the plate 
electrode geometry. 
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It is interesting to note that the charge relaxation time Tq = e^d/as appropriate to thin 
conducting films emerges as the natural unit of time in our analysis. As discussed in Ref. [Q, 
in a thin film the relaxation time is greater than the bulk value e/a by the factor d/ e^s, where 
er is the relative permittivity of the fluid. This is a consequence of the restricted geometry, 
and the fact that the fields lie in the free space outside the film. The wave number of the 
convection pattern observed at onset changes when the film is driven with ac voltages for 
frequencies much larger than l/r^. It would be interesting to modify our analysis to allow 
for time-periodic driving voltages. 

It is often useful to describe patterns near onset with an equation for the slowly varying 
amplitude A of the pattern. For one-dimensional systems which are symmetric under A 
—A, the appropriate amplitude equation is the Ginzburg-Landau equation 

r„f.,A-,\A\'A^(f^. (4.2, 

Here the amplitude A can be taken as the amplitude of the convective velocity field, and 
g governs the nonlinear saturation of the amplitude, and tq are characteristic length 
and time scales introduced earlier. We have previously demonstrated that measurements 
near the onset of convection in smectic films can be described well by Eq. ( ^.2| ) with 
e = (V/Vc)'^ — 1 Ull^l- The onset of convection occurs at a supercritical bifurcation, and the 
e dependence of the flow velocity above onset, the behavior of the amplitude of convection 
near a lateral boundary and the relaxation of the pattern amplitude after sudden changes 
in e are all well described by fits to Eq. ([4.2|) . Eq. ( [4.2|) can also be derived from the full 



electrohydrodynamic equations of motion presented above [|T3 . 

Our analysis gives theoretical results for and tq, as discussed above. Our predicted 
value for correlation length is .^o = 0.285, which is about 20% smaller than the experimental 
value of ^o""^* = 0.36 ± 0.03. This is in fair, but not completely satisfactory, agreement. 
To arrive at both k^^^* and ^o^^*, the experimental measurements were made nondimensional 
by dividing by the measured film width d, which is known to within about 5%. The main 
obstacle to making quantitative comparisons to our predictions for and tq are the 
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poorly known material parameters a and (especially) rj. These appear in the expression 
for the slope V^(s)/s, and are also required for calculating nondimensional times. Realistic 
smectic films are sufficiently viscous that they have values of P ^ f, so that we expect 
the infinite-P limit of tq to apply. The conductivity a for the doped smectic liquid crystal 
used in our experiments has been measured at 1 kHz in a bulk sample |]T^ and at dc in 



an annular film [|l8l. Over this frequency range, it changes by a factor of three. To get 
agreement between Tq^^* and our theoretical value requires a value of a which lies between 
the dc and 1 kHz measurements. Agreement with the Vc{s)/s data |]10|, however, requires 
using a value of r/ a factor of 20 larger than that estimated by extrapolating measurements 



of 7] made in the higher temperature nematic phase [jT9[. This discrepancy may be a result 
of neglecting the effects of air drag on the moving film, which are likely to be important for 
thin, fast-moving films . 



The instability we have described occurs in thin films of fluids which are isotropic in the 
plane of the film. It should also exist for anisotropic fluid films near dc, for example in films 
of smectic C and C* materials. In smectic C materials, the molecules are tilted with respect 
to the layer planes, so the layers are two-dimensional analogs of a nematic fluid. Smectic C* 
materials have an additional broken symmetry which allows a spontaneous electric dipole 
moment in the plane of the layers. Flows in these materials will involve strong orientational 
effects. It should be straightforward to generalize our analysis to the anisotropic case, which 
may lead to interesting new effects. Recently, it has been suggested that electroconvection, 
driven by the analog of the Carr-Helfrich mechanism which operates in negative dielectric 
anisotropy nematics [Q], may occur in smectic C films under ac voltages [^. If materials 
with the right parameters exist, it seems likely that the new instabilities will coexist or 
compete at low frequencies with the instability we have considered here. Something of this 
sort is observed in nematic films [0,0 in which both a "vortex" and a "domain" mode are 
found. 
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V. CONCLUSION 



We have presented a linear stability analysis for the onset of electro convection in a thin 
conducting fluid with two free surfaces. We found the neutral stability curve 71o{k.), along 
with its critical values TZc and Kc, and the correlation length implied by its curvature near 
Kc- The linear growth rate was used to find the characteristic time Tq. This was done for 
two simple electrode configurations. These results were compared with experiments, mainly 
on smectic A films. Several generalizations of this analysis were suggested. 
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TABLES 



TABLE I. Numerical results. 

electrode geometry wires plates 

critical wavenumber, Kc 4.744 4.223 

critical control parameter, Tic 76.77 91.84 

correlation length, 0.2843 0.2975 

time scale, to{V = oo) 0.351 0.372 
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FIGURES 
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FIG. 1. (a) A schematic of the film and electrodes, as used in smectic experiments |^]. The film 
and electrode are shown enormously exaggerated in thickness; in fact s/d ~ 10^^. (b) Schematic 
illustration of the fields inside and outside the film in the small box in part (a), a is a surface 
charge density. 
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FIG. 2. The coordinates used in the analysis, in which the film is treated as a two-dimensional 
sheet. 
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FIG. 3. The surface charge densities for plate (dashed hne) and wire (solid line) electrodes. 
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FIG. 5. The neutral curve for plate (dashed line) and wire (solid line) electrodes. 
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FIG. 6. To as a function of V for wire electrodes. 



